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ABSTRACT 

In the Hubble Ultra Deep Field (HUDF) an abundance of galaxies is seen with a knot at one end 
plus an extended tail, resembling a tadpole. These "tadpole galaxies" appear dynamically unrelaxed 
- presumably in an early merging stage — where tidal interactions likely created the distorted 
knot-plus-tail morphology. Here we systematically select tadpole galaxies from the HUDF and study 
their properties as a function of their photometric redshifts. In a companion HUDF variability study 
presented in this issue, Cohen et al. (2005) revealed a total of 45 variable objects believed to be Active 
Galactic Nuclei (AGN). Here we show that this faint AGN sample has no overlap with the tadpole 
galaxy sample, as predicted by recent theoretical work. The tadpole morphology — combined with the 
lack of overlap with the variable objects — supports the idea that these galaxies are in the process of 
an early-stage merger event, i.e., at a stage that likely precedes the "turn-on" of any AGN component 
and the onset of any point-source variability. We show that the redshift distribution of tadpole galaxies 
follows that of the general field galaxy population, indicating that — if most of the tadpole galaxies 
are indeed dynamically young — the process of galaxy assembly generally kept up with the reservoir 
of available field galaxies as a function of cosmic epoch. These new observational results highlight the 
importance of merger-driven processes throughout cosmic history, and are consistent with a variety 
of theoretical and numerical predictions. 

Subject headings: galaxies: formation — galaxies: mergers — galaxies: active galactic nuclei - 
cosmology 
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1. INTRODUCTION 

The origin of disk galaxies has long been thought to 
occur through the process of dissipational collapse in 
a Cold Dark Matter (CDM) universe (White & Rees, 
1978). Within this paradigm, hierarchical clustering 
(Navarro, Frenk, & White 1997) produces dark matter 
halos in which dissipational collapse of the residual gas 
occurs. The resulting disks retain the kinematic informa- 
tion of their host dark matter potential wells (Blumen- 
thal et al. 1986). Recent numerical simulations have re- 
solved some long-standing discrepancies in the standard 
dissipational collapse scenario by including previously- 
neglected energetic feedback from central supermassive 
black holes during galaxy merging events (e.g. Robertson 
et al. 2005). In particular, they emphasize the relation- 
ship between the central black hole mass and the stellar 
velocity dispersion, which confirms the link between the 
growth of black holes and their host galaxies (di Mat- 
teo, Springel, & Hernquist 2005; Springel, di Matteo, & 
Hernquist 2005ab). These theoretical predictions place 
merger-driven scenarios on the forefront, suggesting that 
galaxy merger activity is a crucial element in a cosmo- 
logical description of the Universe. The present study 
provides observational support for many of these theo- 
retical predictions. 

A large abundance of galaxies in the Hubble Ultra 
Deep Field (HUDF; Beckwith et al. 2005) appear dynam- 
ically unrelaxed, which suggests they must play an im- 
portant role in the overall picture of galaxy evolution. 
In particular, we notice many galaxies with a knot-plus- 
tail morphology. This particular morphology constitutes 
a large, well-defined subset of the irregular and pecu- 
liar objects in the HUDF that is uniquely measurable 



as described in Section 3. The selection of this specific 
morphology also ties closely to the numerical simulations 
described above (di Matteo, Springel, & Hernquist 2005; 
Springel, di Matteo, & Hernquist 2005ab), which predict 
a stage of merger-driven galaxy evolution that closely re- 
sembles these tadpole galaxies in a distinct phase that 
does not yet show AGN activity (as discussed further in 
Section 7). In particular, this morphology appears to 
represent an early stage in the merging of 2 nearly-equal 
mass galaxies. We systematically selected galaxies dis- 
playing this knot-plus-tail morphology from the HUDF; 
a representative sample of tadpoles is shown in Figure ^ 
(details of sample selection are given in Section 3). All 
the selected tadpole galaxies contain the asymmetric, 
pointlike source with a diffuse tail morphology, some with 
multiple knots; all of which we believe are undergoing 
recent interactions. They are mostly linear structures, 
some resembling the "chain" galaxies first reported by 
Cowie, Hu, & Songaila (1995). When more than two 
clumps come together, these objects may be more akin 
to the luminous diffuse objects and clump clusters (Con- 
selice et al. 2004; Elmegreen, Elmegreen, & Sheets 2004; 
Elemegreen, Elmegreen, & Hirst 2004), or other types 
of irregular objects (van den Bergh 2002). Elmegreen 
et al. (2005) visually classify 97 HUDF galaxies (down 
to 10 pixels in size) as "tadpoles" and 126 as "double- 
clump." Some of the galaxies classified by Elmegreen 
et al. as "double-clump" were identified as tadpoles by 
our code, due either to the unresolved nature of one 
clump (which would have been detected as a "tail" in 
our analysis) or to the diffuse nature of one end of the 
object. Since our goal in selecting these tadpoles was to 
sample galaxies that had recently undergone interaction, 
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inclusion of some of these "double-clumps" is warranted. 
One high rcdshift object in our sample has been studied 
in detail by Rhoads et al. (2005). A few objects with mul- 
tiple knots are detected by our selection software, but the 
majority have a single prominant knot with an extended 
tail. 

In this paper, we present the photometric redshift dis- 
tribution of tadpole galaxies, and compare it with the 
rcdshift distribution of the general field galaxy popula- 
tion. This paper is organized as follows. In §2, we de- 
scribe the HUDF data, and in §3 the tadpole sample 
selection. In §4, we discuss why the majority of tadpole 
galaxies are likely not chance alignments, but instead 
mostly dynamically young objects. In §5, we discuss their 
rcdshift distribution, in §6 their relation to galaxy assem- 
bly, and in §7 their possible relation to AGN growth. 

2. HUBBLE ULTRA DEEP FIELD DATA 

The Hubble Ultra Deep Field is a 400 orbit 
survey in four filters carried out using the Ad- 
vanced Camera for Surveys (ACS) aboard the Hub- 
ble Space Telescope (HST) of a single field centered 
on RA(J2000)=03 h 32 m 39?0, Dec(J2000)=-27°47'29'/l. 
The 144-orbit F775W {%') image is deepest, followed by 
F850LP (z'; 144 orbits), F606W (V; 56 orbits), and 
F435W (B; 56 orbits). The HUDF reaches ~1.0 mag 
deeper in B and V and ~1.5 mag deeper in i' than the 
equivalent filters in the Hubble Deep Field (Williams 
et al. 1996). From the ~10,000 objects detected in the 
HUDF (Koekemoer 2004), we will select the sample of 
tadpole galaxies and analyze their properties using the 
i'-band image, because it provides the highest sensitiv- 
ity of the four filters. Yan & Windhorst (2004b) discuss 
how this results in a bias against objects at z > 5.5. This 
bias is small, and only concerns the high redshift tail of 
the redshift distribution. Note, however, that tadpole 
galaxies at z ~ 5.5 do exist (e.g., Rhoads et al.2005). 

3. TADPOLE SAMPLE SELECTION 

The first step in this analysis is to systematically select 
the galaxies that have the characteristic tadpole shape. 
We selected sources in the F775W {%') band to i' AB = 28.0 
mag using SExtractor (Bertin & Arnouts 1996). The 
objects of interest all have a bright "knot" at one end 
with an extended "tail" at the other. SExtractor selects 
objects from an image based on different input parame- 
ters, and adjusting them results in the desired selection 
of sources. The crucial parameter at this stage is DE- 
BLEND_MINCONT, which governs the manner in which 
nearby peaks in flux are considered part of a single ob- 
ject and thus are counted as one source. With the de- 
blending parameter set to a high value, SExtractor will 
separate nearby flux maxima into separate sources. In 
contrast, when the deblcnding is set to a low value, the 
program will count the nearby maxima largely as one 
source. Two different source catalogs are thus generated: 
the highly deblended catalog will contain many point-like 
sources, including the knots of potential tadpole galaxies. 
The catalog with low-dcblcnding will contain extended 
sources, including the tadpole galaxies' tails. The cata- 
logs contain many more sources than the desired ones, 
and the correctly shaped objects must be selected from 
these two initial catalogs. The desired tadpole galaxies 
have a nearly unresolved knot or concentration, and an 



TABLE 1 

Galaxy Selector Input Parameters^ 



Parameter Value 



b/a limit: knots >0.70 

b/a limit: tails <0.43 

Distance to center (in a-axis units) <4 

Angle difference 9 (tail-knot) < 20° 

Total number of tadpoles automatically selected 154 

DEBLENDJV1LNCONT (knots) 0.000005 

DEBLENDJvlINCONT (tails) 0_1 



^Table 1 lists SExtractor parameters used in the tadpole sam- 
ple selection, b & a are semi-major and semi-minor axis lengths 
of the ellipses used in selecting tadpole knots and tails. DE- 
BLENDJVIINCONT (knots/tails) is the SExtractor parameter 
that governs in what manner nearby peaks in flux are considered 
part of a single object and are thus counted as one source. 

extended tail, so these types of sources must be selected 
from the initial catalog, and related spatially such that 
they represent real objects. All of the following proce- 
dures were performed using IDL. 

Both input SExtractor catalogs, described above, con- 
tain the following information about the selected sources: 
x and y pixel locations, length of the semi-major and 
semi-minor axis (a and b) of SExtractor ellipses, and 
the angle (9) of the semi-major axis from north through 
cast. The following input parameters were calibrated us- 
ing a training set of tadpole galaxies that were manually 
selected by visual inspection of portions of the HUDF. 
First, the knots of the tadpole galaxies were selected by 
setting an axis-ratio limit. A "knot" was defined to be 
a source from the highly deblended catalog with an axis 
ratio greater (i.e., rounder) than some critical value (in 
our case, b/a > 0.70). In the same way, the "tails" 
needed to be elongated objects, so a similar procedure 
was performed on the objects from the catalog with low- 
deblending, but with the criterion that their b/a < 0.43. 
The two new lists of correctly shaped objects had to be 
related physically on the image, thus, a new set of objects 
was defined where a knot was within a certain distance of 
the geometrical center of a tail. This distance was taken 
to be < 4a (in semi- major axis units of the tail). We also 
required that the knot be at least > 0.1a from the tail's 
geometrical center, since we are searching for asymmet- 
ric objects, and want to eliminate upfront as many of 
the true edge-on mid-type spiral disks as possible. The 
objects also must have the knot near one end of the tail, 
and this was accomplished by selecting only those tails 
and knots that had a relative angle 9 — measured with 
respect to the semi-major axis of the tail — that was 
< 20°. This step prevented including knots and tails 
that were close together on the image, but not physically 
part of the same galaxy. 

These selection criteria provided a list of the tadpole- 
shaped galaxies. The final number of tadpoles selected 
depended on the selection program's input parameters 
for the limiting axis-ratios, the distance at which the 
knots and tails were considered related, and the angle 
difference between the knot and the semi-major axis of 
the tail. The values of these parameters are given in 
Table 1. With these values, the tadpole galaxy selection 
program detected 154 sources total. This sample was 
then refined as follows. 
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A large majority of the 154 tadpole galaxy candidates 
selected had the characteristic elongated knot-plus-tail 
morphology, although there were some anomalies. In to- 
tal, 14 (9%) obvious mis-detections were visually rejected 
because they were very faint, on the edge of the image, 
or in the outskirts of large face-on spiral galaxies, where 
both knotty and diffuse regions are common and spatially 
close together. A visual examination of the field also pro- 
duced 25 more tadpole galaxies not found by the selector 
program due to the inability of SExtractor to correctly 
separate particular point-like sources within these galax- 
ies. These extra selected tadpole objects visually obeyed 
the morphological criteria that were used to define the 
main sample. Our total final sample thus contains 165 
tadpole galaxies, a subset of which is shown in Figure ^ 
In our final sample, less than 10% of the selected tadpoles 
appear as normal edge-on disk galaxies; the vast major- 
ity have the highly asymmetric morphology. In terms of 
visual vs. automatic selection, we find our sample to be 
about 91% (140/154) reliable and about 86% (154/179) 
complete. The final set of 165 tadpoles galaxies will now 
be studied as a separate class of dynamically unrelaxed 
objects and compared to the general field galaxy popu- 
lation in the HUDF. 

4. WHY TADPOLE GALAXIES ARE NOT CHANCE 
ALIGNMENTS 

In this section we demonstrate that these tadpole 
galaxies are likely not chance alignnments of tails and 
unrelated knots. We first select all elongated diffuse 
structures ("tails") in the HUDF, and then measure the 
angle 8 of the nearest off-centered knot within a radius 
r < 4a (< 2"). Chance alignments of unrelated tails 
and knots would show a random distribution of angles; 
however, Figure El shows that there clearly is an excess 
of knots near \8\ ~ 0°. The excess peak contains 154 
knots, while the average number of knots with \8\ > 10° 
is ~15 per 5°#-bin. Figure [5] thus shows a significant 
overabundance of knots near the end of elongated diffuse 
structures as compared to randomly distributed knots. 
Hence, this physically meaningful result suggests that 
the majority of tadpole galaxies are not just chance align- 
ments of unrelated knots. Instead, we believe they are 
mostly linear structures which are undergoing interac- 
tions. When compared to models of galaxy mergers (di 
Matteo, Springel, & Hernquist 2005; Springel, di Mat- 
teo, & Hernquist 2005; Robertson et al. 2005, Hopkins 
et al. 2005), these objects strongly resemble dynamically 
young objects in the early stages of merging. 

5. THE REDSHIFT DISTRIBUTION OF TADPOLE 
GALAXIES 

To investigate the occurrence of tadpole galaxies 
throughout the history of the universe, we calculate pho- 
tometric redshifts of all HUDF galaxies to i' AB = 28.0 
mag. All photometric redshifts were calculated from 
the HUDF BViz(+JH) photometry using HyperZ (Bol- 
zonella et al. 2000). In order to investigate associated 
redshift errors, we compared our photometric redshifts 
to published spectroscopic redshifts for CDFS 70 ob- 
jects. We find an rms scatter of 0.15 for the fractional 
photometric redshift error (5=(photoz-specz)/(l+specz) 
if all 70 objects are included, and 0.10 when we reject 
a few of the most obvious outliers. This result is fully 



consistent with prior claims of photometric redshift ac- 
curacy in the literature (Lanzetta et al. 1997, Mobashcr 
et al. 2004). The accuracy of our photometric redshift 
estimates depends on the accuracy of the measured mag- 
nitudes in each of the available filters. It also is largely 
independent of the shape of an object (although mag- 
nitude errors for more extended, lower surface bright- 
ness objects tend to be somewhat larger than those for 
more concentrated, higher surface brightness objects of 
the same total magnitude). 

The redshift distribution of all galaxies in the HUDF 
(solid line in Figure |2J is as expected, with the primary 
peak at 0.5 <z< 1.0 and a generally declining tail at 
z~ 4-5. These trends were also seen in the general 
HDF redshift distribution of faint field galaxies (Driver 
et al. 1998). Also apparent is a lack of objects at z~ 1- 
2 due to unavailable UV spectral features crossing the 
BViz(+JH) filters. This occurs because the HUDF 
does not have deep enough F300W or t/-band (ultra- 
violet) data, unlike the situation in the HDF (Williams 
et al. 1996). This redshift bias, however, is the same for 
both the tadpole and the general field galaxy popula- 
tions. In Figure the tadpole galaxy distribution is 
multiplied by a factor of 16 for best comparison with 
that of the field galaxies. Within the available statistics, 
the redshift distribution shape of the tadpole galaxies fol- 
lows that of the general field galaxies quite closely. This 
suggests that if tadpole galaxies are indeed dynamically 
young objects related to early-stage mergers, they may 
occur in the same proportion to the field galaxy popu- 
lation at all redshifts. Tadpole galaxies may therefore 
be good tracers of the galaxy assembly process. The ra- 
tio of the two redshift distributions N(z) was calculated 
as well, and the resulting percentage of tadpole galaxies 
is plotted in Figure 01 as a function of redshift together 
with the statistical errors. Overall, the percentage of tad- 
pole galaxies is roughly constant at ^6% with redshift to 
within the statistical errors for the redshift range probed 
in our study (0.1 <z< 4.5). 

6. TADPOLE GALAXIES AS TRACERS OF GALAXY 
ASSEMBLY 

The fact that about 6% of all field galaxies are seen in 
the tadpole stage is a measurement with potentially im- 
portant consequences. In light of simulations by Springel 
et al. (2005) that predict a tadpole-like stage ~0.7 Gyr 
after a major merger begins, we suggest that this particu- 
lar tadpole morphology represents an early-merger stage 
of two galaxies with comparable mass. If this 6% indi- 
cates the fraction of time that an average galaxy in the 
HUDF spends in an early-merger stage during its life- 
time, and if most of these low-luminosity objects started 
forming the bulk of their stars at the end of the reion- 
ization epoch at z~ 6 — 7 (e.g. Yan & Windhorst 2004a, 
2004b), then each galaxy would spend about 6% of 12.9 
Gyr (i.e. 0.8 Gyr) since z~ 7 in a distinctly recognizable 
merger or tadpole stage. At the median redshift at which 
the tadpoles are seen (z me( j ~ 1.6; see FigureEJ, each ob- 
ject is then seen at an age of about 4 Gyr if born at z~7. 
Each tadpole is ~ 1" (or ~ 8 kpc) across (Figure^]), and 
given the fluxes measured, each clump in a tadpole has 
roughly M ~ 1O 8 -1O 9 M0 in stars (see, e.g., Papovich, 
Dickinson, & Ferguson 2001 who estimated stellar masses 
of lyman break galaxies). For these rough estimates of 
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Fig. 1. — F775W («') band mosaic of a subset of the tadpole galaxy sample in the HUDF. Stamps retain the orientation of the HUDF; 
north is toward the top of the page and east is to the left. Index numbers are displayed in the upper-left corner of the stamps; photometric 
redshifts are given in upper-right corners. Parentheses indicate photometric redshifts with errors >1 based on HyperZ calculations. Stamps 
are 3 arcsec on a side. A table of coordinates for the entire tadpole sample is given in Section 3; the entire sample of 165 tadpoles appears 
in color in the online supplement. The vast majority of our tadpole sample contains the distinctive knot-plus-tail morphology, while sample 
contamination by normal (non-interacting) edge-on disk galaxies is less than 10%. 
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Fig. 2. — Distribution of (the absolute value of) angles (8) of all 
off-centered knots found within a radius r < 4a (< 2") from the 
center of each elongated diffuse structure in the HUDF, showing a 
clear excess of knots near |6*| ~ 0°. 

their physical parameters, the freefall timescale for each 
tadpole is roughly t < (fewx 10 7 )-10 8 years, or ~6% 
of the galaxy lifetime at that redshift. Hence, if every 
galaxy is seen in a tadpole stage for ~0.8 Gyr of its life- 
time, then it may have undergone ^10-30 mergers dur- 
ing its lifetime. During the early stage of each merger, 
it would be temporarily seen as a tadpole. More com- 
plex mergers may lead to irregular/peculiar and train- 
wreck type objects and the luminous diffuse objects or 
clump clusters, which are among the type of objects that 
dominate the galaxy counts at faint magnitudes (Driver 
et al. 1998). In this paper, we limit the sample selection 
to two clumps passing by each other, which we believe 
leads to the more uniquely classifiable tadpole morphol- 



FlG. 3. — Photometric redshift distribution of galaxies in the 
HUDF. The solid black histogram shows the redshift distribution of 
all HUDF field galaxies to i' AB = 28.0, while the dashed histogram 
shows the redshift distribution of the tadpole galaxies. The latter 
was multiplied by 16 X for best comparison of its shape with the 
redshift distribution of the field galaxies. 

ogy. Given that tadpoles only trace a certain type and 
stage of merging galaxies, the above statistics are likely 
a lower limit on the number of all mergers. In summary, 
each galaxy seen today may have had of order one- to 
two- dozen mergers since most of its Population II stars 
were born at z~7, and given the small masses and short 
merger timescales involved, would then be seen as tad- 
pole galaxies for about 6% of their life-time. 

Figure 0] suggests that tadpole galaxies — if indeed 
dynamically young objects — appear to occur in the 
same proportion to the field galaxy population at all 
redshifts probed in this study. Tadpole galaxies may 
therefore be good tracers of the process of galaxy as- 
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Fig. 4. — Percentage of total galaxies that arc tadpoles is plotted 
as a function of photometric redshift. Within the statistical errors, 
on average about 6% of all galaxies are seen as tadpoles at all 
redshifts. 

scmbly. This implies that the process of galaxy assembly 
- as traced by tadpole galaxies — keeps up with the 
reservoir of available galaxies as a function of redshift 
for 0.1 < z < 4.5. Our result is in excellent agreement 
with the predictions of Robertson et al. (2005) that de- 
scribe a merger-driven scenario to build up disk galaxies, 
and is consistent with Rhoads et al. (2005) who conclude 
that their z=5.4 galaxy (tadpole 165 in Table 2) is 
strongly indicative of a galaxy in assembly. 



7. 



RELATION BETWEEN TADPOLE GALAXIES AND 
ACTIVE GALACTIC NUCLEI 



In a companion paper in this issue, Cohen et al. (2005) 
present a study of the variable objects in the HUDF 
which have a point source component with a measur- 
ably variable flux on timescales of 0.4-3.5 months, which 
roughly corresponds to 0.5-5.5 weeks in the rest-frame 
at the median redshift of the sample (z me d ~ 1.5 — 2). 
In particular, they found 45 plausibly variable objects 
among 4644 galaxies to i' AB = 28.0 mag. They argue that 
these objects are most likely variable because they host 
weak AGN. Sometimes these AGN are in the galaxy cen- 
ter, but often they occur off-center in a dynamically unre- 
laxed system. This prompts the question: What fraction 
of tadpole galaxies contains a variable weak AGN in its 
knots? This is a critical issue, because it is widely be- 
lieved that the process of merging in galaxies can also 
disturb the inner accretion disk around the supermas- 
sive black hole (SMBH) and switch on the AGN. Among 
our 165 tadpole galaxies, none coincide with the sam- 
ple of 45 variable objects or with the x-ray sources in the 
Chandra Deep Field South (Alexander, D.M. et al. 2005). 
Recent state-of-the-art hydrodynamical models (di Mat- 
teo, Springel, & Hernquist 2005; Springel, di Matteo, & 



Hernquist 2005; Hopkins et al. 2005) suggest that during 
(major) mergers, the black hole accretion rate peaks con- 
siderably after the merger started, and after the star for- 
mation rate (SFR) has peaked. Specifically, their models 
suggest that, for massive galaxies, a tadpole stage is seen 
typically about 0.7 Gyr after the merger started, but 
^0.9 Gyr before the SMBH accretes most of its mass, 
which is when the galaxy displays strong AGN activ- 
ity. Since the lifetimes of QSO's and radio-galaxies are 
known to be < (fewx 10 7 )-10 8 years (Martini, P. 2004, 
Grazian et al. 2004), these hydrodynamical models thus 
imply that the AGN stage is expected to occur consid- 
erably after (i.e., >1 Gyr) the early-merger event dur- 
ing which the galaxy is seen in the tadpole stage. The 
observed lack of overlap between the tadpole galaxies 
and the AGN sample in the HUDF provide direct ob- 
servational support for this prediction. Recent studies 
by Grogin et al. (2005) find asymmetry (A) values that 
are similar between AGN and non-AGN samples; this re- 
sult is consistent with our study in light of the theoretical 
models mentioned above, which indicate that AGN activ- 
ity is seen only well after the merger has taken place and 
the galaxy has settled into a more dynamically relaxed 
state. In addition, Hopkins et al. (2005) have quantified 
the timescales that quasars will be visible during merging 
events, noting that for a large fraction of the accretion 
time, the quasar is heavily obscured. In particular, their 
simulations show that during an early merging phase — 
our observed tadpole phase — the intrinsic quasar lu- 
minosity peaks, but is completely obscured. Only after 
feedback from the central quasar clears out the gas will 
the object become visible as an AGN. This should be 
observable by the Spitzer Space Telescope in the mid- 
infrared (IR) as a correspondingly larger fraction of IR- 
selected obscured faint QSO's. 

In conclusion, tadpole galaxies are a class of easily- 
identifiable, dynamically young objects that exist 
throughout the history of the Universe and are good 
tracers of galaxy assembly. They provide strong obser- 
vational support for the validity of recent numerical sim- 
ulations, and highlight the importance of mergers to the 
process of galaxy assembly and AGN growth. 
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Global Properties of Tadpole Galaxies 
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39 
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07 
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45 
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35 
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16 
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46 
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35 
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4t 



NOTE: parentheses indicate estimated aperture magnitudes for visually selected objects 
^rcdshift from Rhoads et al. 2005 

t Photometric redshifts computed from HyperZ (Bolzonella et al. 2000) 



